Classical, semiclassical, and quantum investigations of the 4-sphere scattering system 
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A genuinely three-dimensional system, viz. the hyperbolic 4-sphere scattering system, is inves- 
tigated with classical, semiclassical, and quantum mechanical methods at various center-to-center 
separations of the spheres. The efficiency and scaling properties of the computations are discussed 
by comparisons to the two-dimensional 3-disk system. While in systems with few degrees of freedom 
modern quantum calculations are, in general, numerically more efficient than semiclassical methods, 
this situation can be reversed with increasing dimension of the problem. For the 4-sphere system 
with large separations between the spheres, we demonstrate the superiority of semiclassical versus 
quantum calculations, i.e., semiclassical resonances can easily be obtained even in energy regions 
which are unattainable with the currently available quantum techniques. The 4-sphere system with 
touching spheres is a challenging problem for both quantum and semiclassical techniques. Ifere, 
semiclassical resonances are obtained via harmonic inversion of a cross-correlated periodic orbit 
signal. 

PACS numbers: 05.45.-a, 03.65.Sq 
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I. INTRODUCTION 

The breakthrough for the semiclassical quantization of 
chaotic systems was the development of periodic orbit 
theory 0j Gutzwiller's trace formula the density 

of states is expressed as an infinite sum over all isolated 
periodic orbits of the classical system. Although the pe- 
riodic orbit theory is in principle valid for systems with 
an arbitrary number of degrees of freedom, applications 
have, for practical reasons, so far mostly been restricted 
to two-dimensional systems. The main difficulties are, 
firstly, the numerical periodic orbit search, which be- 
comes more difficult in multidimensional systems, and, 
secondly, the fact that the semiclassical trace formula 
usually does not converge. The convergence problems 
can be solved, e.g^ wit_h cycle-expansion 0, 0, Q or har- 



monic inversion [a Pi techniques, and both methods 
have been successfully applied to the 3-disk billiard as a 
prototype model of a two-dimensional hyperbolic scatter- 
ing system. Practical applications of periodic orbit the- 
ory to three-dimensional systems are very rare. For the 
three-dimensional Sinai billiard extensive quantum com- 
putations have been performed and the quantum spectra 
have been analyzed in terms of classical periodic orbits 
0,^3- However, no semiclassical eigenstates have been 
calculated from the set of periodic orbits. Semiclassical 
resonances have been obtained for the three-dimensional 
2- and 3-sphere scattering systems but for these sys- 
tems all periodic orbits lie in a one- or two-dimensional 
subspace. 

In this paper will investigate the scattering of a par- 
ticle on four equal spheres centered at the corners of 
a regular tetrahedron. Classical and semiclassical as 
well as quantum mechanical methods will be applied 
to the 4-sphere system at various center-to-center sep- 
arations of the spheres. The 4-sphere system can be 
regarded as the simplest extension of the 3-disk repel- 



lor to three-dimensional space with a set of genuinely 
three-dimensional periodic orbits. Chaotic properties of 
the 4-sphere system have been verified experimentally 
by the observation of fractal structures via optical light- 
scattering on the spheres 0, 0| . 

When solving two- and three-dimensional systems with 
both quantum and semiclassical methods it is interesting 
to study the scaling properties of the quantization meth- 
ods with the number of degrees of freedom, and to com- 
pare the efficiency of the various algorithms. The numer- 
ical effort for the quantization of nonintcgrable systems 
usually increases strongly with the number of degrees of 
freedom, and therefore efficient quantization techniques 
are highly desirable. A large variety of quantum me- 
chanical and semiclassical methods have been developed. 
The direct solution of Schrodinger's equation is possi- 
ble, e.g., by time-dependent wave packet expansions or 
numerical diagonalization of the Hamiltonian in a com- 
plete basis set. Exact quantum mechanical calculations 
usually require storage of multidimensional wave func- 
tions and a computational effort that grows exponen- 
tially with the number of coupled degrees of freedom. 
These methods are therefore feasible for systems with 
relatively few degrees of freedom. As an alternative to 
exact quantum calculations, approximate, e.g. semiclas- 
sical, methods can be used. Gutzwiller's trace formula 
can be applied to systems with an arbitrary number of 
degrees of freedom, however, the number of periodic or- 
bits and the numerical effort needed to find them usually 
increases very rapidly with increasing dimension of the 
phase space. As a matter of fact, Gutzwiller's periodic 
orbit theory has been applied almost exclusively to sys- 
tems with two de gree s of freedom, viz. the anisotropic 
Kepler problem [lLll4| , the hydrogen atom in a magnetic 



field [15j, and two-dimensional bilhards 3, 4, 16, 17j. 
For these systems direct quantum mechanical computa- 
tions are usually more powerful and efficient than the 
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semiclassical calculation of spectra by means of periodic 
orbit theory. The 4-sphere system is an example where 
semiclassical methods turn out to be superior to direct 
quantum mechanical computations |l8l |. i.e., semiclassi- 
cal resonances can easily be obtained even in energy re- 
gions which are unattainable with the presently known 
quantum techniques. 

The paper is organized as follows. In Sec. we in- 
vestigate the classical dynamics of the 4-sphere system. 
The symbolic code is introduced and its symmetry re- 
duction by means of the tetrahedra group, Td, is dis- 
cussed. The periodic orbits are found in a systematic 
way by an efficient numerical periodic orbit search, and 
the pruning of orbits at small separations of the spheres 
is analyzed. In Sec. IIIII we introduce the semiclassical 
techniques for periodic orbit quantization, viz. the cycle- 
expansion method, the harmonic inversion method, and 
the extension of harmonic inversion to cross-correlated 
periodic orbit signals. In Sec. llVl we present the method 
applied for the exact quantum mechanical calculation of 
the resonances. In Sec. |3 we show the results for the 
semiclassical and quantum resonances at various separa- 
tions of the spheres. The results are discussed with spe- 
cial emphasis on the comparison of the efficiency of the 
various methods. Some concluding remarks are given in 
Sec. Ell 



II. CLASSICAL DYNAMICS: THE PERIODIC 
ORBITS OF THE 4-SPHERE SYSTEM 

The 4-sphere system is a genuinely three-dimensional 
billiard where the systematic periodic orbit search is a 
nontrivial task. In this section we first develop the sym- 
bolic dynamics of orbits and the symmetry reduction us- 
ing the tetrahedra group, and then discuss the numerical 
periodic orbit search and the calculation of the periodic 
orbit parameters. 



A. Symbolic code and symmetry reduction 

The 4-sphere system discussed here consists of four 
equal spheres with radius a centered at the corners of 
a regular tetrahedron. We choose a = 1 in what fol- 
lows. The system is then solely determined by the center- 
to-center-separation R. The 4-sphere system with large 
center-to-center-separation R ^ 2a and with touching 
spheres {R — 2a) are shown in Fig.n(a) and (b), respec- 
tively. 

In full coordinate space each orbit can be described by 
the infinite sequence of spheres where the orbit is scat- 
tered. By labehng the spheres as {A, B, C, D}, it is possi- 
ble to code a periodic orbit as the infinite cycles of a lim- 
ited length string consisting of the sphere labels which we 
call here the itinerary code of the orbit. For a given string 
length, all combinations of the letters {A, B,C, D} cor- 
respond to a physical orbit, with the exception that two 
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FIG. 1: The 4-sphere system consists of four equal spheres 
centered at the corners of a regular tetrahedron, (a) Large 
center-to-center-separation R S> 2a. (b) Touching 4-sphere 
system with R = 2a. 



consecutive letters in the itinerary code cannot be iden- 
tical and, for short center-to-center-separation R > 2a, 
some orbits may be excluded by pruning (see Sec lIICfl . 
Several itinerary code strings may represent the same 
periodic orbit or a similar orbit obtained by a symme- 
try operation, i.e., rotation or reflection. For example, 
the itinerary codes ABC and BCA correspond to the 
same periodic orbit by cyclic permutations, and the or- 
bits ABC, ACD, ABD, and BCD can be mapped onto 
each other by rotations. 

By using the symmetry properties Td of the tetrahe- 
dron the system can be reduced to its fundamental do- 
main. The symmetry reduced orbits can be described by 
a ternary alphabet of symbols '0', '1', and '2', which are 
the three fundamental orbits, i.e., the symmetry reduc- 
tions of the shortest orbits scattered between two, three, 
and four spheres, respectively. Therefore, we shall use the 
symbol '0' for returning back to the previous sphere after 
one reflection, symbol '1' for the reflection to the other 
third sphere in the same reflection plane of the orbit, and 
symbol '2' for the reflection to the other forth sphere out 
of the reflection plane of the orbit. The reflection plane is 
defined by the centers of the first three different spheres 
toward back in the history of the itinerary code of the 
orbit. The primitive periodic orbits of cycle length rip in 
the fundamental domain are now given by those periodic 
sequences of rip symbols '0', '1', and '2' which are free 
of subcycles (e.g. the code 0101 with subcycle 01 is not 
primitive, we will neglect the line indicating periodicity 
in the following). The periodic orbits do not change by 
cyclic permutations of the code. We will choose the code 
word with the lowest numerical value as the representa- 
tive (e.g. 0112 instead of 1120). With these rules every 
symmetry reduced periodic orbit of the 4-sphere system 
is uniquely described by a symbolic code. However, at 
small separations of the spheres some physical orbits are 
pruned as discussed below in Sec. Ill CI 

From the {0,1,2} code of the symmetry reduced or- 
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TABLE I: Symbolic code p of the symmetry reduced periodic 
orbits with cycle lengths < 3 and the itinerary codes p of 
the orbits in full coordinate space. The column hp gives the 
symmetry type of the orbits. 



p 


P 


hp 





AB 




1 


ABC 




2 


ABDC 


S4, 


01 


ABAC 




02 


ABADAC 


C3 


12 


ABO JJBAIJC 


64 


001 


ABABCBCAC 


C3 


002 


ABABDBDCDCAC 


S4 


Oil 


ABACBC 


CTd 


012 


ABACDC 


C2 


021 


ABADBDCBC 


Cs 


022 


ABADCDBABCDC 


S4 


112 


ABCADC 




122 


ABCDACBDC 


C3 



bits the {A, B, C, D} itinerary code can be obtained as 
follows. We choose the plane spanned by the spheres 
{A, B, C) as the initial reflection plane and start the jour- 
ney with the sequence AB. Then the rules given above 
are applied for the symbols 0, 1, and 2 to guide the or- 
bit to the subsequent spheres. Note that symbolic codes 
which contain only the symbols '0' and '1' lie in the two- 
dimensional (A, B, C)-plane, i.e., they correspond to the 
set of orbits with a binary symbolic code, which has been 
well-established for the 3-disk 0, 3 and 3-sphere system 
pT]. Orbits including the symbol '2' are genuinely three- 
dimensional orbits. In Table we present the symbolic 
codes of all periodic orbits up to cycle lengths n^, = 3 of 
the symmetry reduced code. Note that no subcycles and 
cyclic permutations exist on the list. In the second col- 
umn, the corresponding itinerary codes p of the symbolic 
codes of column 1 are given, which have been obtained 
by following the rules explained above. The last column 
in Table ^ shows the symmetry classes of the orbits. The 
Td group has 1 e, 3 C2, 8 C3, 6 S4, and 6 ad, in total 
24 different symmetry elements. Each orbit (except the 
one represented by 0) can be assigned by one and only 
one of the symmetry elements {e, (Xd, C2, C3, 6*4} of the 
group Td- Note that periodic orbits in the fundamen- 
tal domain, and thus their symmetry reduced symbolic 
codes, are two-, three-, or four-times shorter than the or- 
bits (and the itinerary codes) in the full coordinate space 
when they belong to the symmetry class {ad, C2}, C3, or 
54, respectively. The symbolic length of orbits belonging 
to symmetry class e, i.e., the identity is unchanged under 
symmetry reduction. 



B. Numerical periodic orbit search 

Each trajectory of the 4-sphere system is completely 
determined by the reflection points on the surfaces of the 
spheres, which on each sphere can be described by two 
spherical coordinates 9 and 0. For a given itinerary code 
arbitrarily chosen reflection points on the spheres con- 
nected by straight lines in the correct order result in a 
periodic but not necessarily a physical orbit. The true 
physical orbit, for which the incident and reflection angle 
at each reflection point must coincide, can be obtained by 
direct application of Hamilton's principle, i.e., the orbital 
length, which is proportional to the classical action, be- 
comes a minimum when the reflection points are varied. 
The length function of an orbit with a total number of 
N reflection points depends on the 2N variables {9i,(j)i} 
with i — 1, . . . ,N. Numerically, the minimizing of the 
length 

L{0i,(j)i,02,h,---,ON,M (1) 

can be achieved by applying the well established quasi- 
Newton method |19i |. which is implemented, e.g., in the 
NAG-library The required gradient of the length 

function, VL, has been derived analytically. 

As mentioned above the length of periodic orbits in 
full coordinate space can be two, three, or four times 
the length of the corresponding symmetry reduced or- 
bit in the fundamental domain (see Table QJ. As the re- 
quired computational effort for the quasi-Newton method 
increases rapidly with the dimensionality of the problem, 
it is desirable to exploit the symmetry properties of the 
tetrahedra group and to directly search for the periodic 
orbits in the fundamental domain. To this end for a sym- 
metry reduced orbit with cycle length Up the reflection 
point on the sphere Up + l is associated with the reflection 
point on the flrst sphere by an appropriate transforma- 
tion, i.e., one of the 24 possible symmetry transforma- 
tions of the tetrahedra group, Td- The length minimiza- 
tion is now applied to the trajectory segments between 
the flrst sphere and sphere Up + 1, i.e., the dimensional- 
ity of the length minimization of periodic orbits in the 
fundamental domain is reduced to 2np for all primitive 
orbits with cycle length Up. 

Once a periodic orbit has been found its orbital pa- 
rameters required for semiclassical periodic orbit quanti- 
zation can be calculated. The most important ones are 
the monodromy matrix and the Maslov index of the orbit. 
The Maslov index increases by 2 at each reflection on a 
hard sphere, i.e., /ipo = 2np for an orbit with cycle length 
Up. The calculation of the monodromy matrix Mpo for 
the periodic orbits of three-dimensional billiards has been 
investigated in Refs. 0, |^. Mpo is a symplectic (4 x 4) 
matrix with eigenvalues Ai, 1/Ai, A2, and I/A2. For the 
hyperbolic 4-sphere system Ai and A2 are either both 
real or the orbits are loxodromic, i.e., the eigenvalues of 
Mpo are a quadruple {A, 1/A, A*, 1/A*} with A being a 
complex number. For the 4-sphere system with radius 
a = 1 and center-to-center separation R = 6 the orbital 
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TABLE 11: Parameters of the symmetry reduced primitive 
periodic orbits p with cycle length < 3 of the 4-sphere 
system with radius a = 1 and center-to-center separation R = 
6. 



p 


hp 


L 


Re Ai 


Im Ai 


Re A2 


Im A2 







4.000000 


9.89898 


0.00000 


9.89898 


0.00000 


1 


C3 


4.267949 


-11.7715 


0.00000 


9.28460 


0.00000 


2 


S4 


4.296322 


-4.52562 


9.49950 


-4.52562 


-9.49950 


01 




8.316529 


-124.095 


0.00000 


88.4166 


0.00000 


02 


C3 


8.320300 


-37.1479 


98.0419 


-37.1479 


-98.0419 


12 




8.567170 


117.644 


0.00000 


-102.992 


0.00000 


001 


C3 


12.321747 


-1240.54 


0.00000 


868.915 


0.00000 


002 


Si 


12.322138 


-353.853 


976.176 


-353.853 


-976.176 


Oil 


Od 


12.580808 


1449.55 


0.00000 


824.981 


0.00000 


012 


C2 


12.617350 


1192.83 


0.00000 


-1020.66 


0.00000 


021 


C3 


12.584068 


1201.43 


0.00000 


-996.800 


0.00000 


022 


54 


12.619948 


-755.582 


804.976 


-755.582 


-804.976 


112 




12.835715 


-496.339 


1038.46 


-496.339 


-1038.46 


122 


C3 


12.863793 


-1100.56 


0.00000 


1219.28 


0.00000 



lengths and stability parameters for all primitive periodic 
orbits with cycle length rip < 3 are presented in Table ITTI 
For that sphere separation (i? = 6a) we have calculated 
the complete set of primitive periodic orbits with symbol 
lengths Up < 14, numbering 533830 orbits in total. For 
sphere separation R — 2.5 a we also calculated all prim- 
itive periodic orbits with symbol lengths Tip < 14, and 
in addition all orbits with symbol lengths Up < 22 and 
physical lengths L < 12, which allows for the construc- 
tion of a periodic orbit signal with length L^ax = 12 used 
for the semiclassical quantization in Sec. lVBl 

C. Pruning of orbits 

For center-to-center separations R > 2.0482 a between 
the spheres there is a one to one correspondence between 
the symbolic codes and the primitive periodic orbits. 
However, when the separation is reduced below that value 
some orbits become unphysical, i.e., the symbolic dynam- 
ics is pruned. The pruning of orbits has been investigated 
in detail for the 3-disk scattering system p3.l2^. For the 
4-sphere system the mechanism is similar: As illustrated 
in Fig. 121 an orbital segment may (a) pass through one of 
the spheres, or (b) a reflection may occur inside one of the 
spheres. For a periodic orbit search at small separation 
between the spheres all orbits obtained numerically by 
minimizing the length must be checked whether pruning 
occurs or not. For touching spheres (i? = 2a) all pruned 
orbits with symbol lengths Up < 7 and their pruning 
types (a) or (b) are presented in Table Hill Pruning exists 
for orbits with symbol lengths Up > 5, i.e., the symbolic 
dynamics is complete only for Up < 4. Furthermore, peri- 
odic orbits with long heads of '0' symbols in the code can 




(a) (b) 



FIG. 2: Sketch of the two types of pruning occurring in the 
4-sphere system at small separation between the spheres; (a) 
An orbital segment passes through one of the spheres, (b) A 
reflection occurs inside a sphere. 



TABLE III: All pruned orbits with cycle length Up < 7 and 
their pruning types (a) or (b) (see Fig. I^J of the 4-sphere 
system with touching spheres, R — 2a. 



symbolic code 


pruning type 


00021 


a 


000011 


a 


000021 


a 


000002 


b 


0000001 


b 


0000011 


a 


0000021 


a 


0000002 


b 



have accumulation points at finite values of the physical 
length L. It is impossible to find all orbits beyond the 
first accumulation point. We have searched for all peri- 
odic orbits of the touching 4-sphere system with physical 
lengths L < Lmax = 3.6, symbol lengths rip < 60, and 
with the total number of '1' and '2' symbols in the sym- 
bolic code restricted to ni + n2 < 10, resulting in about 
2.8 million primitive periodic orbits. 

The semiclassical quantization by harmonic inversion 
of a cross-correlated periodic orbit signal (see Sec. lIII (T)) 
requires the knowledge of the expectation values of vari- 
ous linearly independent classical observables A along the 
periodic orbits 24, 25] . We have chosen the observables 
Ai — and A2 = L'^, i.e., we have calculated the aver- 
aged squared distance and squared angular momentum 
of the periodic orbits of the touching 4-sphere system. 



III. SEMICLASSICAL PERIODIC ORBIT 
THEORY 

We now wish to calculate semiclassically the reso- 
nances of the 4-sphere scattering system by application 
of periodic orbit theory. Gutzwiller's trace formula 
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TABLE IV: Weight factors Wp for the periodic orbit sum 
(Eq. (2J) of the 4-sphere system with symmetries of the tetra- 
hedra group, Td- 





e 




C2 


Si 




^1 


1 


1 


1 


1 


1 


A2 


1 


1 


1 


-1 


-1 


E 


2 


-1 


2 








Ti 


3 





-1 


1 


-1 


T2 


3 





-1 


-1 


1 



expresses the quantum mechanical response function 

in terms of the periodic orbits of the underlying classical 
system, i.e.. 



{E)/h 



(3) 



po 



where go{E) is a smooth function of the energy and 
Apo(-E) and Spo{E) are the periodic orbit amplitudes (in- 
cluding phase information given by the Maslov indices) 
and classical actions, respectively. For billiards the classi- 
cal action depends linearly on the length of the trajectory 
and the wave number k = \/2ME/h with M being the 
particle mass. For the three-dimensional 4-sphere system 
the periodic orbit sum as a function of the wave number 
k reads 



5W = EE 



'LpC 



ikrLr} 



r=l 



(4) 

where Hp is the cycle length, Lp the physical length, Xp i 
are the eigenvalues of the monodromy matrix, and r is 
the repetition number of the primitive periodic orbit p. 
The weight factors Wp result from the symmetry decom- 
position of the system |26| and depend on the chosen 
irreducible subspace of the spectrum and the symmetry 
of the periodic orbits. For the tetrahedra group, Td, the 
values of the weight factors Wp are given in Table IIVI 
In the following we will concentrate on the subspace Ai , 
where the weight factors of all orbits are Wp — 1. 

The semiclassical resonances of the 4-sphere system are 
given by the poles of the function g{k) . However, it is well 
known that the periodic orbit sum Q does not converge 
in those regions where the physical poles are located, and 
special techniques must be applied to obtain an analytical 
continuation of the periodic orbit sum Q . For the 3-disk 
system with large center-to-center separation R = 6a the 
cycle-expansion method [^I^ l27j and harmonic inversion 
techniques (6|, |^ have proven to be powerful approaches 
for overcoming the convergence problems of the periodic 



orbit sum, and both methods can also be successfully ap- 
plied to the 4-sphere system. However, when pruning of 
orbits sets in at small separations, and in particular in the 
case of touching disks or spheres, the situation is much 
more difficult and subtle, since the direct application of 
the cycle-expansion method fails. The two-dimensional 
closed 3-disk billiard is a bound system, where a few semi- 
classical eigenenergies have been obtained in Ref. us- 
ing the cycle-expansion in combination with a functional 
equation. This method is not valid for open systems and 
cannot be extended to the 4-sphere system which remains 
open even in the case of touching spheres |0, 0| . Nev- 
ertheless, the harmonic inversion of cross-correlated pe- 
riodic orbit signals 0j IHl has been successfully applied 
to the closed 3-disk system '28', '2^ and this method will 
also serve as a powerful tool for the three-dimensional 4- 
sphere system. We will now introduce the quantization 
methods. Applications to the 4-sphere system and com- 
parisons with quantum mechanical results will be pre- 
sented in Sec. IVl 



A. The cycle-expansion method 

The periodic orbit sum in Gutzwiller's trace formula 
does usually not converge in the energy regions of phys- 
ical interest. However, for some systems, e.g., the 3-disk 
scattering billiard, semiclassical energies or resonances 
can be obtained with the help of the cycle-expansion 
method ^ 4, 5j. If the periodic orbits can be associated 
to a symbolic dynamics the Gutzwiller-Voros zeta func- 
tion [l|,|30| can be expanded according to increasing cycle 
length of the orbits. In this expansion the contributions 
of long periodic orbits may be approximately shadowed 
by the combined contributions of shorter orbits. In this 
case the cycle-expansion can converge rapidly. 

For billiards the Gutzwiller-Voros zeta function can be 
written as 



Zovik; 2) = exp <! - ^ ^ - 

p r—l 



1 (—zY^fe"''''"^ 



det(Mj; 



(5) 



with an additional parameter z which must be set to 
z = 1. The cycle-expansion is achieved by taking z as a 
book-keeping variable and expanding Eq. ((Sj as a trun- 
cated power series in z. The semiclassical resonances 
are obtained as the zeros (in the variable k) of the cycle- 
expanded zeta function (O with again z — 1. In our com- 
putations for the 4-sphere system we use cycle-expansions 
up to order n,jiax — 12. 



B. Semiclassical quantization by harmonic 
inversion 

An alternative method for semiclassical quantization is 
based on the observation that the extraction of eigenval- 
ues from Gutzwiller's trace formula can be reformulated 
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as a signal processing task 0, S ■ The harmonic inver- 
sion method is briefly explained as follows. The Fourier 
transform of the function g{k) in Eq. Q yields the semi- 
classical signal 



det(Mr - 1)1 



-J{L - rLp) , (6) 



as a sum of 5 functions. The central idea of semiclassical 
quantization by harmonic inversion is to adjust the semi- 
classical signal C^^{L) with finite length L < Lmax to its 
quantum mechanical analogue 



1 

2^ 



+ 00 



^ k 

n 



dn 



ikL 



dk 



(7) 



where the amplitudes d„ and the semiclassical eigenval- 
ues kn are free adjustable co mplex p arameters. This is 
achieved by signal processing [sj, jsj of the semiclassical 
signal C^'^(L). Numerical recipes for extracting the pa- 
rameters {dm kn} by harmonic inversion of the 5 function 
signal © are given in P.l33l|. 



C. Harmonic inversion of cross-correlated periodic 
orbit signals 

The method of semiclassical quantization by harmonic 
inversion of cross-correlated periodic orbit signals is a 
generalization of the quantization scheme presented in 
Sec. IIII Bl The idea is to use the classical average values 
of a set of linearly independent classical observables to 
construct a cross-correlated signal, whose informational 
content is significantly increased as compared to the one- 
dimensional signal, and therefore should lead to semiclas- 
sical spectra with improved resolution. 

The numerical tools for the harmonic inversion of 
cross-correlated periodic orbit signals have already been 
well established (251 , and therefore we only briefly review 
the basic ideas and refer the reader to the literature for 
details. For simplicity but without loss of generality, we 
focus on billiard systems, where the shape of the orbits 
is independent of the energy E, and the classical action 
of orbits reads S — hkL, with k the wave number and L 
the physical length. The starting point is to introduce a 
weighted response function in terms of k 



9aa' {k) 



(8) 



where fc„ is the eigenvalue of the wave number of eigen- 
state \n) and 



ban {n\Aa\n) 



(9) 



are the diagonal matrix elements of a chosen set of N 
linearly independent operators Aa, a = 1,2^ . . . , N. The 



Fourier transform of (jSj) yields the NxN cross-correlated 
signal 



-t-oo 



Cac'iL) = — / 5«„,(fc)e-"=^dfc 

ZTT , 



(10) 



A semiclassical approximation to the cross-correlated sig- 
nal H10|l has been derived in [23 . l83 |. The cross-correlated 
periodic orbit signal reads 



^Q:,p '^o'.p( l) ^ Lip 



r=i ,/|det(M;-l) 



d{L - rLp) , 



(11) 

where r is the repetition number counting the traversals 
of the primitive orbit, and Mp is the monodromy matrix 
of the primitive periodic orbit. The weight factors Ua^p 
are classical averages over the periodic orbits 



A„(q(L),p(L))dL, (12) 



JO 



with Aa (q, p) the Wigner transform of the operator Aa ■ 
Semiclassical approximations to the eigenvalues fc„ and 
eventually also to the diagonal matrix elements (n|AQ,|n) 
are now obtained by adjusting the semiclassical cross- 
correlated periodic orbit signal 1)11(1 to the functional 
form of the quantum signal H1U|) . The numerical tool 
for this procedure is an extension of the harmonic inver- 
sion method to the signal processing of cross-correlation 
functions [s^ [s^. The advantage of using the cross- 
correlation approach is based on the realization that the 
total amount of independent information contained in 
the NxN signal is 7V(A^ + 1) multiphed by the length 
of the signal, while the total number of unknowns (here 
ban and fc„) is {N + 1) times the total number of poles 
kn- Therefore the informational content of the NxN 
signal per unknown parameter is increased (as compared 
to the one-dimensional signal) by roughly a factor of TV, 
and the cross-correlation approach should lead to a sig- 
nificant improvement of the resolution. 



IV. QUANTUM CALCULATIONS 

Schrodinger's equation for the 3-disk or the 4-sphere 
system is a free wave equation in two or three dimen- 
sions, [A -I- /c^]^'(k) = 0, with Dirichlet boundary con- 
ditions, i.e., ^(k) = on the surface of the disks or 
spheres, respectively. Although the problem looks simple 
the solution is a nontrivial task and, most importantly, 
the numerical effort increases extremely rapidly with the 
dimension of the system. 

For the 3-disk system the exact quantum resonances 
can be obtained as roots of the equation |22i, l37ij 



detM(fc) 



3-disk 



0, 



(13) 
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with m and m! nonzero integer numbers which can be 
truncated by an upper angular momentum TOmax ^ 
1.5 ka 27]. With matrices 'M.{k)'^^}^ of dimension up to 
^ (400 X 400), Eq. 1^3) allows for the cfhcicnt numerical 
calculation of resonances in the region < Refca < 250. 
The matrix elements of M(A;)^^';'^ in Eq. lO can be 
written analytically in terms of Bessel and Hankel func- 
tions. Explicit expressions are given in 27]. The quan- 
tum resonances are obtained by a numerical root search 
in the complex fc-plane |l9l ]. 

Similarly, exact quantum resonances of the three- 
dimensional 4-sphere scattering system can be obtained 
as roots of the equation 



detM(fc)f-Pt'^ = 



(14) 



with < Z, < Z,„ax and m, m' — 0, 3, 6, 9, ... , Zmax for 
the subspace Ai and A2. An explicit expression for the 
matrix elements of '^{k)^'^^^'^ has been derived ^ and 
reads 



jiika) 
h[!\kay 



Po)hfikR)il5) 



with 



C{l,m,l',m'J;0o,l3o) 



I' 

E 

M=-l' 



{2l + l){2l' + l){2l + l) 



I V I 




(-1)"51,™-m(^o,0) 
±^:-™-m(^o,0) 



m ~ M M m 



(16) 



where the it-signs refer to the subspace Ai and A2^ re- 
spectively. The angles 6*0 and /3o in (fTB|l are obtained 
from 



Note that 5m=o should read go = l/\/2 instead of in 
Eq. (38) of [ol. Similar as for the 3-disk system the an- 
gular momentum (quantum numbers I and V in Eq. (fTHl l 
can be truncated at Z^ax ^ 1.5 fca to achieve convergence 
of the calculation. 

It is important to note that the computation of the 
quantum mechanical resonances of the three-dimensional 
4-sphere scattering system becomes much more expen- 
sive than for the two-dimensional 3-disk system. First of 
all, the calculation of each matrix element lV[(fc)^^^^'^[,° 
in (|15|l requires the summation over quantum numbers 
I and (via (|16|l ) M . To accelerate the calculation of the 
matrix H15|l at various values of k we have calculated and 
stored the values of C{l,m,l' ,m' ,1;9q, Pq) in (|16|) sepa- 
rately. Eq. Hlt)|l does not depend on the sphere separa- 
tion R, and therefore the stored C-values can be used in 
calculations of spectra with arbitrary R. However, the 
calculation of the matrix elements in (|15|l still requires 
the summation over I. 

The second problem of solving Eq. 114|) is the scaling 
of the dimension of the matrix '^{k)f^^^^? , which is an 
N X N matrix with 



N : 



^ (^max 




-3), 



i.e., N scales as N ^ k"^ for the 4-sphere system, as com- 
pared to N ^ k for the 3-disk system, Eq. H13(l . For 
example, in the region ka ~ 200 the required matrix di- 
mension is > 300 for the 3-disk system, as compared to 
N > 15000 for the 4-sphere system. For the 4-sphere sys- 
tem with center-to-center separations R — 6a, R — 2.5 a, 
and the touching spheres R — 2a we have computed 
the quantum resonances in the region < Re ka < 60 
by solving Eq. H14|l with matrices of dimension up to 
(1751 X 1751). The results wih be presented in Sec. |V| 
With currently available computer technology it is im- 
possible to significantly extend the quantum calculations 
for the 4-sphere system to the region Re fca ^ 60 using 
Eq. (|14|) . The efficiency of the semiclassical and quantum 
methods for the 4-sphere system will be compared and 
discussed in Sec. IV Dl 



cos(6'o) 
cos(/3o) 



2 

Ve 
1 

3 ' 



sin (6*0) 
sin(/3o) 



1 

|V2, 



and the di„^,(d) are the matrix elements of finite rota- 
tions 1381. 



The large brackets in H16|l refer to 3j-symbols 
the values of Qm are defined as 

r l/\/2 form = 
9m ^ ■( 1 for m = 3, 6, 9, . . . , Z 

1 otherwise. 



and 



V. RESULTS AND DISCUSSION 

We will now present and discuss the results of our semi- 
classical and quantum computations for the 4-sphere sys- 
tem with large sphere separation R = 6a, intermediate 
separation R = 2.5 a, and touching spheres, R = 2a. In 
Sec. IV Dl we will compare and discuss the efficiency of the 
various quantization methods. 



A. Sphere separation R — 6a 

The quantum mechanical and semiclassical Ai- 
resonances of the 4-sphere system with radius a = 1 
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FIG. 3: Ai-resonances in the complex fc-plane of the 4-sphere 
system with radius a = 1 and center-to-center separation 
R = 6. Squares: Quantum computations. Crosses and 
plus symbols: Semiclassical resonances obtained by cycle- 
expansion and harmonic inversion methods, respectively. 



and center-to-center separation i? = 6 are presented in 
Fig. 121 The quantum resonances marked by the squares 
have been obtained by solving Eq. (|14(l with matrices 
M{k)f-^lf^? of dimension up to (1134 x 1134), which is 
sufficient only to obtain converged results in the region 
Re fc < 50 (see Fig. By contrast, the semiclassical 
resonances can easily be obtained in a much larger region, 
e.g.. Re fc < 250 shown in Fig. The crosses mark the 
zeros of the cycle-expanded Gutzwiller-Voros zeta func- 
tion |SJ| . The cycle-expansion has been truncated at cycle 
length riniax = 7, which means that a total set of just 508 
primitive periodic orbits are included in the calculation. 
The plus symbols mark the semiclassical resonances ob- 
tained by harmonic inversion of the periodic orbit signal 
© with signal length Lmax = 60 constructed from the 
set of 533830 primitive periodic orbits with cycle lengths 
Up < 14. 

In the region Re fc < 50 (Fig. |3Ji) the quantum and 
semiclassical resonances agree very well, with a few ex- 
ceptions. The first few quantum resonances in the up- 
permost resonance band are narrower, i.e., closer to the 



real axis than the corresponding semiclassical resonances. 
A similar discrepancy between quantum and semiclassi- 
cal resonances has already been observed in the 3-disk 
system 0, Furthermore, in the region Re fc < 15 

and Im k < —0.5 several quantum resonances have been 
found (see the squares in Fig.|3i), which seem not to have 
any semiclassical analogue. These resonances are related 
to the diffraction of waves at the spheres, and its semi- 
classical description requires an extension of Gutzwiller's 
trace formula and the inclusion of diffractive periodic or- 
bits 1131 • The semiclassical resonances obtained by 
either harmonic inversion or the cycle-expansion method 
(the plus symbols and crosses in Fig. respectively) 
are generally in perfect agreement, except for the very 
broad resonances that lie deep in the complex plane, i.e., 
in the region Im k < —0.8. 



B. Sphere separation R = 2.5 a 

The semiclassical quantization becomes more and 
more demanding with decreasing separation between the 
spheres. The reason is that the shadowing of longer orbits 
by combinations of shorter orbits in the cycle-expanded 
Gutzwiller-Voros zeta function becomes less accurate and 
the construction of the periodic orbit signal of length 
L < imax used for the harmonic inversion method re- 
quires more and more periodic orbit data. However, 
both semiclassical quantization techniques, i.e., cycle- 
expansion and harmonic inversion can still be success- 
fully applied at significantly reduced separation between 
the spheres. 

As an example of an intermediate sphere separation we 
discuss the case R = 2.5 a, where the spheres are rather 
close, however, the symbolic dynamics of the periodic 
orbits is still complete, i.e., no orbits are pruned. The 
graphical comparison of the quantum mechanical and 
semiclassical resonances in the region < Re fca < 100 
is given in Fig. 0] The semiclassical resonances shown 
as plus symbols have been obtained by harmonic inver- 
sion of a periodic orbit signal of length L,nax = 12. The 
signal has been constructed using all primitive periodic 
orbits with symbol lengths Up < 14 and parts of the or- 
bits with symbol lengths 15 < Up < 22, in total a set 
of about 4.6 million orbits. The crosses in Fig. 0] mark 
the semiclassical resonances obtained by 12th order cycle- 
expansion using the complete set of 69706 primitive pe- 
riodic orbits with symbol lengths Up < 12. The exact 
quantum resonances have been obtained in the region 
< Re fc < 60 by solving Eq. H14|) with matrix dimen- 
sions up to (1751 X 1751). 

In the region Re fc < 60 the resonances obtained by 
the two semiclassical methods are in excellent agreement 
except for the imaginary parts of some resonances very 
deep down in the complex plane. In this region the 
semiclassical resonances agree well with the exact quan- 
tum mechanical resonances, the deviations are due to the 
semiclassical approximation, i.e., the first-order h expan- 
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FIG. 4: Ai-resonances in the complex fc-plane of the 4-sphere 
system with radius a = 1 and center-to-center separation 
7? = 2.5. Squares: Quantum computations. Crosses and 
plus symbols: Semiclassical resonances obtained by cycle- 
expansion and harmonic inversion methods, respectively. 



sion in the semiclassical trace formula. As in the case 
R — 6a (Sec. IV A|l some quantum resonances in the re- 
gion Re fc < 10 are related to the diffraction of waves 
at the spheres and do not have a semiclassical analogue 
vifithout the ap propriate extension of the periodic orbit 
theory [s^ ^ . At Re fc > 60 the agreement between 
resonances obtained semiclassically via cycle-expansion 
and harmonic inversion becomes less perfect, especially 
for some broad resonances with Im fc < —1.1. Unfor- 
tunately, no quantum results are currently available for 
Re fc > 60 to judge the quality and accuracy of the semi- 
classical computations in that region. 



C. Four touching spheres {R = 2a) 

The semiclassical quantization of the 4-sphere system 
becomes even more difficult when the spheres are fur- 
ther moved together and the symbolic dynamics becomes 
pruned (see Sec. Ill C|) . In particular, the case of touching 
spheres with i? = 2a is a real challenge for the follow- 
ing reason. For touching spheres the symbolic dynam- 
ics is pruned in a similar way as in the 3-disk problem 
[23|. The closed 3-disk billiard is a bound system, and 
some eigenenergies have been extracted by either combin- 
ing the cycle-expansion method with a functional equa- 
tion (14] or by the harmonic inversion method 2^ i2ft | . 
However, contrary to the closed 3-disk system the four 
touching spheres do not form a bound system, which 
means that the method of Ref. 14^ combining the cycle- 
expansion method with a functional equation cannot be 
applied, and thus the touching 4-sphere system cannot be 
quantized with the help of the cycle-expansion method. 
Nevertheless, we will now demonstrate that the harmonic 
inversion method applied to a cross-correlated periodic 



FIG. 5: Ai-resonances in the complex fc-plane of the touch- 
ing 4-sphere system with radius a = 1 and center-to- 
center separation R = 2. Squares: Quantum computations. 
Crosses: Semiclassical resonances obtained by harmonic in- 
version without using cross-correlation. Plus symbols: Semi- 
classical resonances obtained by the harmonic inversion of a 
(3 X 3) cross-correlated periodic orbits signal using the oper- 
ators 1 (identity), L^, and r^. 



orbit signal can reveal at least some of the low-lying semi- 
classical resonances. 

For the construction of the periodic orbit signal we 
have calculated about 2.8 million orbits of the touching 4- 
sphere system with lengths L < L^^i,^ = 3.6. (Note that 
the signal is incomplete as discussed in Sec. Ill CI ) For 
the application of the cross-correlation technique we use 
the operators Ai = 1 (the identity), the squared angular 
momentum A2 = Lp' , and the squared distance from the 
origin, = . Because the signal is incomplete and 
rather short the results of the harmonic inversion are 
less perfectly converged than for the 4-sphere system with 
larger sphere separation, i.e., the amplitudes d„ in Eq. 
may deviate from the ideal values dn — ^ for true physical 
resonances and d„ = for spurious resonances which 
must be omitted. As a criterion to accept resonances we 
have chosen the condition |c?„ — 1| < 0.5. 

The results of our semiclassical and quantum computa- 
tions for the four touching spheres are presented in Fig.|Sl 
The crosses mark the semiclassical resonances obtained 
by harmonic inversion of the one-dimensional periodic or- 
bit signal. The low number of crosses indicates that the 
convergence of the one-dimensional signal is not very sat- 
isfactory. The plus symbols show the resonances obtained 
by harmonic inversion of the (3 x 3) cross-correlated pe- 
riodic orbit signal using the operators Ai = 1, = L^, 
and A3 = r^. With the cross-correlation technique the 
convergence properties have been significantly improved 
compared to the analysis of the one-dimensional signal. 
The real parts of the semiclassical resonances agree well 
with the real parts of the exact quantum mechanical res- 
onances marked by the squares in Fig. |S1 The agreement 
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between the imaginary parts is, however, less perfect. 
Some quantum resonances in Fig. [S] do not have a semi- 
classical counterpart. Those resonances with Re /c < 10 
are probably related to the diffraction of waves at the 
spheres as discussed above, i.e., they cannot be explained 
without extensions of the semiclassical theories applied in 
this paper. 



D. Efficiency of tlie semiclassical and quantum 
algorithms 

For the 4-sphere system as an example of a physical 
system with three degrees of freedom we now wish to 
discuss and compare the efficiency of the semiclassical 
and quantum computations. As mentioned in the intro- 
duction (Sec. nj the efficiency of quantum computations 
usually decreases rapidly with the number of degrees of 
freedom of the physical system. It is an interesting and 
important question whether semiclassical methods can 
beat the efficiency of quantum computations with in- 
creasing dimension of the problem. Although there is 
not much hope and evidence that this is generally true, 
because of the exponential proliferation of periodic or- 
bits in chaotic systems, it can be true for certain specific 
systems. An example for the superiority of semiclassi- 
cal over quantum mechanical calculations is the 4-sphere 
system with large sphere separation, e.g. R = 6a, where 
semiclassical resonances can easily be obtained even in 
energy regions which are out of reach for the presently 
known quantum techniques |l8l |. To understand this it 
is instructive to study the expense and scaling properties 
of the quantum and classical computations for the 3-disk 
and 4-sphere system. 

As explained in Sec. IIVI exact quantum resonances of 
the 3-disk and 4-sphere systems can be obtained as roots 
of Eqs. (|13|l and (|14|l . respectively, with angular quan- 
tum numbers truncated at Imax ^ 1.5 ka. The calcula- 
tion of the matrix elements in H14|l is much 
more expensive than for the matrix elements M(fc)^j^/j'5'^ 
in H13(l . However, the serious problem of solving Eq. H14(l 
is the scaling of the dimension of the matrix M(fc)^^']!'^", 
which is an iV X matrix with = (Zmax+2)(Zmax+3)/6, 
i.e., A^ scales as A^ fc'^ for the 4-sphere system, as 
compared to N ^ k for the 3-disk system, Eq. (|13|l . 
For example, in the region ka ~ 200 the required ma- 
trix dimension is A^ > 300 for the 3-disk as compared 
to A^ ^ 15000 for the 4-sphere system. With currently 
available computer technology it is, therefore, impossible 
to significantly extend the quantum calculations for the 
4-sphere system to the region ka 60 using Eqs. ifCT 
I16|l . Note that the cost of the quantum computations 
does not depend on the separation R between the disks 
or spheres. 

The expense of the semiclassical quantization is ba- 
sically given by the required number of periodic orbits 
which, in chaotic systems, increases exponentially with 



the symbolic or physical length of the orbits. For the 3- 
disk system the number of symmetry reduced primitive 
periodic orbits with symbol length Up is given approxi- 
mately by Af ~ 2"p / Up whereas it scales as A^ ^ 3"p /up 
for the 4-sphere system. Contrary to the quantum com- 
putations the numerical expense for the semiclassical 
quantization, i.e., the required number of orbits depends 
on the separation R between the disks or spheres. For 
large separation R = 6a the cycle-expansion method is 
most efficient for the calculation of a large number of 
resonances. The reason is that the assumption of the 
cycle-expansion that the contributions of longer peri- 
odic orbits in the expansion of the Gutzwiller-Voros zeta 
function Q are shadowed by pseudo-orbits composed of 
shorter periodic orbits is very well fulfilled. The har- 
monic inversion method also allows for the calculation of 
a large number of resonances, but requires a larger input 
set of periodic orbits. While for the two-dimensional 3- 
disk system the semiclassical and quantum computations 
are very efficient, the semiclassical methods are superior 
to the quantum techniques for the three-dimensional 4- 
sphere system. The semiclassical calculations can easily 
be extended to the region Re ka > 60 where no quantum 
results are available because of the unfavorable scaling 
of the dimension of the matrix M;m,i'm' in Eq. (|14|l . Of 
course, a more efficient quantum method for the 4-sphcre 
system than that of Rcf. IL] may in principle exist. How- 
ever, to the best of our knowledge no such method has 
been proposed in the literature to date. The 4-sphere 
system therefore is an example of a three-dimensional 
system where semiclassical methods are presently supe- 
rior to exact quantum calculations. 

At reduced separation R = 2.5 a between the disks 
or spheres the semiclassical quantization requires an in- 
creased set of periodic orbits to achieve convergence of 
the cycle-expansion or harmonic inversion analysis. How- 
ever, for the 4-sphere system the semiclassical methods 
are still superior to the exact quantum computations, 
i.e., semiclassical resonances can be obtained in regions 
which are unattainable with the quantum methods as can 
be seen in Fig. 0] 

The situation is different for touching spheres, R — 2a, 
which is a challenging system not only for the quantum 
but also for the semiclassical computations. The con- 
struction of a long periodic orbit signal is impossible be- 
cause orbits with increasing sequences of consecutive '0' 
symbols in the code lead to accumulation points in the 

tsical length similar as for the closed 3-disk system 
|2^ . The semiclassical calculations for the touching 
spheres are therefore at least about the same or even 
more expensive than the quantum computations. 



VI. CONCLUSIONS 

In summary, we have investigated an open sys- 
tem with three degrees of freedom, viz. the 4-sphere 
scattering problem with various sphere separations by 
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means of classical, semiclassical, and quantum mechan- 
ical methods. The classical system has genuinely three- 
dimensional periodic orbits. In the symmetry reduced 
fundamental domain, they can be associated to a ternary 
symbolic alphabet, which allows for a systematic periodic 
orbit search. For large separations between the spheres 
{R > 2.5 a) semiclassical resonances have been obtained 
by application of the cycle-expansion technique and the 
harmonic inversion method. For touching spheres (i? = 
2a), the symbolic dynamics is pruned and the cycle- 
expansion does not converge, however, some semiclas- 
sical resonances can be revealed by harmonic inversion 
of a cross-correlated periodic orbit signal. 

Exact quantum mechanical resonances have also been 
calculated, however, the quantum computations for the 
three-dimensional 4-sphere system are much more expen- 
sive than for the two-dimensional analogue, viz. the 3- 
disk scattering problem. Therefore, the quantum com- 
putations had to be restricted to the region with rela- 
tively low wave numbers, i.e., Refca < 60. By analyzing 
the scaling properties of both the quantum and semiclas- 
sical calculations we have demonstrated the superiority 
of semiclassical methods over quantum computations at 
least for large sphere separations, i.e., semiclassical reso- 
nances can easily be obtained in energy regions which at 
present are unattainable with the established quantum 
method. These results may encourage the investigation 
of other systems with three or more degrees of freedom 



with the goal of developing powerful semiclassical tech- 
niques, which are competitive with or even superior to 
quantum computations for a large variety of systems. 

In those regions where exact quantum results for the 4- 
sphere system are lacking an assessment of the accuracy 
of the semiclassical resonances is presently impossible. 
Higher-order h corrections have been calculated for two- 
dimensional billiard systems [4ll li^ . l43l |. however, the 
extension of the theory to three-dimensional systems is a 
nontrivial task for future work. 

Those quantum resonances which are related to diffrac- 
tion of waves at the spheres have not yet been explained 
semiclassically. For the 3-disk system diffractive reso- 
nances have been obtained with an extended periodic 
orbit theory by including the contributions of creeping 
orbits [33, • It will be interesting to generahze these 
ideas to the genuinely three-dimensional 4-sphere system. 
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